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Poincare covariant Faddeev equations for the nucleon and A are solved to illustrate that an 
internally consistent description in terms of confined-quark and nonpointlike confined-diquark- 
correlations can be obtained. vrA'^-loop induced self-energy corrections to the nucleon's mass are 
analysed and shown to be independent of whether a pseudoscalar or pseudovector coupling is used. 
Phenomenological constraints suggest that this self-energy correction reduces the nucleon's mass 
by up to several hundred MeV. That effect does not qualitatively alter the picture, suggested by 
the Faddeev equation, that baryons are quark-diquark composites. However, neglecting the vr-loops 
leads to a quantitative overestimate of the nucleon's axial- vector diquark component. 
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I. INTRODUCTION 

Contemporary experimental facilities employ large mo- 
mentum transfer reactions to probe the structure of 
hadrons and thereby attempt to elucidate the role played 
by quarks and gluons in building them. Since the pro- 
ton is a readily accessible target its properties have been 
studied most extensively Q. Hence an understanding of 
a large fraction of the available data requires a Poincare 
covariant theoretical description of the nucleon. 

At its simplest the nucleon is a nonperturbative three- 
body bound-state problem, an exact solution of which 
is difficult to obtain even if the interactions are known. 
Hitherto, therefore, phenomenological mean-field models 
have been widely employed to describe nucleon structure; 
e.g., soliton models ||, ^ and constituent-quark mod- 
els 1^, 1^, . These models are most naturally applied to 
processes involving small momentum transfer (q^ < Af^, 
M is the nucleon mass) and, as commonly formulated, 
their applicability may be extended to processes involv- 
ing larger momentum transfer by working in the Breit 
frame |8||. Alternatively, one could define an equivalent, 
Galilean invariant Hamiltonian and reinterpret that as 
the Poincare invariant mass operator for a quantum me- 
chanical theory Q but this path is less well travelled. 

Another approach is to describe the nucleon via a 
Poincare covariant Faddeev equation. That, too, requires 
an assumption about the interaction between quarks. An 
analysis Q of the Global Colour Model |^, |l| sug- 
gests that the nucleon can be viewed as a quark-diquark 
composite. Pursuing that picture yields a Faddeev 
equation, in which two quarks are always correlated as a 
colour-antitriplet diquark quasiparticle (because ladder- 
like gluon exchange is attractive in the 3c quark-quark 
scattering channel) and binding in the nucleon is effected 
by the iterated exchange of roles between the dormant 



and diquark-participant quarks. 

A first numerical study of this Faddeev equation for 
the nucleon was reported in Ref . , and following that 
there have been numerous more extensive analyses; e_^ 
Refs. [^6[ |T^. In particular, the formulation of Ref. 
employs confined quarks, and confined, pointlike-scalar 
and -axial-vector diquark correlations, to obtain a spec- 
trum of octet and decuplet baryons in which the rms- 
deviation between the calculated mass and experiment is 
only 2%. The model also reproduces nucleon form fac- 
tors over a large range of momentum transfer [p^ , and its 
descriptive success in that application is typical of such 
Poincare covariant treatments; e.g., Refs. ||l9|, |2^, p2| . 

However, these successes might themselves indicate a 
flaw in the application of the Faddeev equation to the nu- 
cleon. For example, in the context of spectroscopy, stud- 
ies using the Cloudy Bag Model (CBM) |||| indicate that 
the dressed-nucleon's mass receives a negative contribu- 
tion of as much as 300-400 MeV from pion self-energy 
corrections; i.e., SM+ = -300 to -400 MeV §, ||. Fur- 
thermore, a perturbative study, using the Faddeev equa- 
tion, of the mass shift induced by pointlike-7r exchange 
between the quark and diquark constituents of the nu- 
cleon obtains 5M+ = -150 to -300 MeV H]. Uname- 
liorated these mutually consistent results would much di- 
minish the value of the 2% spectroscopic accuracy ob- 
tained using only quark and diquark degrees of freedom. 

It is thus apparent that the size and qualitative impact 
of the pionic contribution to the nucleon's mass may pro- 
vide material constraints on the development of a real- 
istic quark-diquark picture of the nucleon, and its inter- 
pretation and application. Our article is an exploration 
of this possibility and we aim to clarify the model de- 
pendent aspects. We emphasise, in addition, that chiral 
corrections to baryon magnetic moments and charge radii 
are also important |23] , and their model- independent fea- 
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tures furnish additional constraints on any quark model, 
including those based on the Faddeev equation, thereby 
guiding their improvement. We note, too, that lattice- 
QCD studies of baryon masses, especially as a function 
of the current-quark mass [ p6| , also provide information 
that can guide these considerations; e.g., a recent lattice- 
QCD exploration of the connection between TV and A 
masses is consistent with the pion self-energies described 
above 

In Sec. II we recapitulate on the Faddeev equation 
and its solution for the N and A in a simple model. 



masses (isospin symmetry is assumed): 



Section II] discusses mo del- independent aspects of the 
Dyson-Schwinger equation (DSE) that describes the 
pionic correction to the iV's self-energy and therein we 
also present exemplary estimates for the magnitude of 
the effect. Section ^ is an epilogue. 
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The mass ordering is characteristic and model- 
independent (cf. Refs. |3^, lattice-QCD estimates |3^ 
and studies of the spin-flavour dependence of parton dis- 
tributions [^), and indicates that a study of the N and 
A must retain at least the scalar and pseudovector [uu)- 
and (M(i)-correlations if it is to be accurate. (Of course, 
the spin-3/2 A is inaccessible unless pseudovector corre- 
lations are retained.) 



Model for the Nucleon 



II. FADDEEV EQUATION 

The properties of light pseudoscalar and vector mesons 
are well described by a renormalisation-group-improved 
rainbow-ladder truncation of QCD's DSEs ||, |^, ll], 
and the study of baryons via the solution of a Poincare 
covariant Faddeev equation is a desirable extension of 
the approach. The derivation of a Faddeev equation for 
the bound state contribution to the three quark scat- 
tering kernel is possible because the same kernel that 
describes mesons so well is also strongly attractive for 
quark-quark scattering in the colour-antitriplet channel 
(see Sec. [I A 2). And it is a simple consequence of the 
Clebsch- Gordon series for quarks in the fundamental rep- 
resentation of SUc{i)- 



> 3c = (3c e 6c) «) 3c = Ic 



10c, (1) 



that any two quarks in a colour singlet bound state must 
constitute a relative colour antitriplet. This supports a 
truncation of the three-body problem wherein the inter- 
actions between two selected quarks are added to yield 
a quark-quark scattering matrix, which is then approxi- 
mated as a sum over all possible diquark pseudoparticle 
terms: Dirac-scalar -I- -pseudovector +[...]- essentially a 
separable two-body interaction |^ . A Faddeev equation 
follows, which describes the three-body bound-state as a 
composite of a dressed-quark and nonpointlike diquark 
with an iterated exchange of roles between the dormant 
and diquark-participant quarks. The bound-state is rep- 
resented by a Faddeev amplitude: 



(2) 



where the subscript identifies the dormant quark and, 
e.g., ^1^2 are obtained from "if^ by a correlated, cyclic 
permutation of all the quark labels. 

The Faddeev equation is simplified further by retain- 
ing only the lightest diquark correlations in the repre- 
sentation of the quark-quark scattering matrix. A sim- 
ple, Goldstone-theorem-preserving, rainbow-ladder DSE- 
model p3] yields the following diquark pseudoparticle 



To provide a concrete illustration and make our pre- 
sentation self-contained we consider a simple model |2^ 
wherein the nucleon is a sum of scalar and pseudovector 
diquark correlations: 



(4) 



with (pi, ai,Ti) the momentum, spin and isospin labels of 
the quarks constituting the nucleon, and P = pi + P2 + 
P3 the nucleon's total momentum. The scalar diquark 
component in Eq. (^ is 

^-f (K,a„T,) = 

[T'\lp[ur,K)]ZZ A°'(A1 [Sii;PHP)]Z , 

(5) 

where [^: the spinor satisfies 

{i"fP + M) u{P) = = u{P) {ij-P + M) , (6) 

with M the mass obtained in solving the Faddeev equa- 
tion, and is also a spinor in isospin space with ip+ = 
col(l,0) for the proton and = col(G, 1) for the neu- 
tron; if = pi -f p2 =: P{12}, P[12] = Pi ~ P2, t ■■= 
(~P{12} + 2p3)/3; A° {K) is a pseudoparticle propagator 
for the scalar diquark formed from quarks 1 and 2, and 
r° is a Bethe-Salpeter-like amplitude describing their 
relative momentum correlation; and 5, a 4 x 4 Dirac 
matrix, describes the relative quark-diquark momentum 
correlation. (5, r° and A*^ are discussed below.) The 
pseudovector component is 

* {p^,ai,n) = 



[evl (-P[i2];if)]?:4 K.iK) [A^,{t-P)u{P)] 



where the symmetric isospin-triplet matrices are 



T3 

as ' 



= -(r° 



(7) 



(8) 
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with {T^)ij = 5ij and t^'^ the usual Pauh matrices, and 
the other elements in Eq. (Q) are obvious generalisations 
of those in Eq. (^) . 

The colour antisymmetry of is implicit in T"^ , with 
the Levi-Civita tensor, ecic2C3 , expressed via the antisym- 
metric Gell-Mann matrices; i.e., defining 

{H^ =iX\H^ = -iX'^,H^ = iX^}, (9) 

then e,,,,,3 = (i/^^),,,,. (See Eqs. (||), (11).) 

The Faddeev equation satisfied by '^^ yields a set of 
coupled equations for the matrix valued functions S, A^: 



S{k;P)u{P) 
Aiik;P)u{P) 



M{k,l- P) 



S{£;P)u{P) 
Aii£;P)uiP) 



(10) 



where one factor of "2" appears because '^3 is coupled 
symmetrically to 5'i and ^'2, and we have evaluated the 
necessary colour contraction: {H'^)bc{H'')ch' = —2Sbb'- 
The kernel in Eq. M) is 



with 



Mik,e] P) 



Mm 
(Mio) 



{Moi)i 



(11) 



^qq/'^] ^qq) S {iq 



{i, - k,j2- -k,,) S{Q , (12) 



The general forms of 5 and A^^, the Bethe-Salpeter- 
like amplitudes that describe the momentum-space cor- 
relation between the quark and diquark in the nucleon, 
are discussed at length in Ref. [1^, wherein a detailed 
analysis of the Faddeev equation's solution is presented. 
Requiring that S be an eigenfunction of A+(P), Eq. (|A^ ), 
entails 



Si£;P)^ hii;P)h 



j^{tj-£-£-PIu) f2{£;P), 

(18) 



where {Iu)rs ~ ^rs, = ^1, and, in the nucleon rest 
frame, /i^2 describe, respectively the upper, lower compo- 
nent of the bound-state nucleon's spinor. Requiring the 
same of A^^ reduces to only six (from an original twelve) 
the number of independent Dirac amplitudes required to 
specify it completely. However, we simplify this by re- 
taining only those two amplitudes that survive in the 
non-relativistic limit: 



AM- P) ^ al(£; P) 757m + a\{£- P) 757 ' «m - (-^ 



1. 



(19) 



1,2. 



Assuming isospin symmetry then: 0} = < ^ 

The Faddeev equation for the nucleon is Eq. ( p^ ) with 
the kernel, A^, given by Eqs. (|ll|)-([T5|): to complete its 
definition we must specify the dressed-quark propagator, 
the diquark Bethe-Salpeter amplitudes and the diquark 
propagators. 



where: £q = £ + P/3, fc, = fc + P/3, £qq = -£+'2Pl'i, kqq 

—k + 2P/3; S is the propagator of the dormant dressed- 
quark constituent of the nucleon (Sec. IIAl); and 



{Moi)i = t^' T]^ (fc, - £qj2- £qq) S^{£qq - fc,) 



xf°"(£. 



l5(£,)Ai;(^,,), (13) 



(Xio); = rO^(A:, - iqj2;£qq) S^{£qq - k^) 



(ylq fcgg/2, ■ 



-fc,,)^(£,)AO {£qq), (14) 



{M,^i,=t^Tl {kq-£qq/2-£qq)S^{£. 



xt'fi"(£. 



kqq/'^'i kqq) ^i^q) ^ 



i^q 



pv V-qq) 



(15) 



In Eqs. (10)- (Em) it is implicit that u{P) is a normalised 



average of ip± so that, e.g., the equation for the proton 
is obtained by projection on the left with ip^j^. To clarify 
this, by illustration, we note that Eq. ( |l^ generates an 
isospin coupling between u{P)^^ on the l.h.s. of Eq. dTT 
and, on the r.h.s.. 



y^2Atu{P)^- -Alu{P) 



(16) 



This is merely the Clebsch-Gordon coupling of isospin- 
1® isospin-i to total isospin-i and means that the scalar 
diquark amplitude in the proton, (ud)Q+ u, is coupled to 
itself and the linear combination: 



y/2 {uu)i+ d — {ud)i+ u . 



(17) 



1. Dressed-quarks 

The general form of the dressed-quark propagator is 

S{p) = -i-fp<Jv{p^) + <Js{p''), (20) 
= [z7-pA(p2) + B(p2)]-i. (21) 

It can be obtained by solving the QCD gap equation; i.e., 
the DSE for the dressed-quark self energy, and the many 
such studies 
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4C|| yield the model-independent re- 
sult that the wave function rcnormalisation and dressed- 
quark mass: 



Z{p') = lM(p2) , M{p') = B{p')/A{p^) 



(22) 



respectively, exhibit significant momentum dependence 
for ^ 1 GeV^, which is nonperturbative in origin. This 
behaviour was recently observed in lattice-QCD simula- 
tions and Refs. [^2|, ^ provide quantitative com- 
parisons between those results and a modern DSE model. 
The infrared enhancement of M{p'^) is an essential conse- 
quence of dynamical chiral symmetry breaking (DCSB) 
and is the origin of the constituent-quark mass. With 
increasing the mass function evolves to reproduce the 
asymptotic behaviour familiar from perturbative anal- 
yses, and that behaviour is unambiguously evident for 
p2 > lOGeV^ PI. 
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While numerical solutions of the quark DSE are readily 
obtained, the utility of an algebraic form for S{p) is self- 
evident. An efficacious parametrisation of S{p), which 
exhibits the features described above, has been used ex- 
tensively in studies of meson properties ^ and we 
use it herein. It is expressed via 

asix) = 2mT{2{x + fh^)) 

+T{hix) Tihx) [bo + b2T{ex)] , (23) 

^v{x) = -i^ [l-.F(2(x-Hm2))] , (24) 



X + 



with X = p^/X^ 



T{x) = [1 — exp(— x)]/x. 



erg (a;) ~ Acrs(p^) and crvix) = A^cry(p^). The mass- 
scale, A = 0.566 GeV, and parameter values 



m 



0.00897 0.131 2.90 0.603 0.185 



(25) 



were fixed in a least-squares fit to light-meson observ- 
ables P4|. The dimensionless u = d current-quark mass 
in Eq. (p5[) corresponds to 



m = 5.1 MeV. 



(26) 



(e = 10""' in Eq. (p3| ) acts only to decouple the large- 
and intermediate-p^domains.) 

The parametrisation expresses DCSB, giving a Eu- 
clidean constituent-quark mass 



Mi 



0.33 GeV, 



(27) 



defined as the solution of p^ — M^(p^), whose 
magnitude is typical of that employed in constituent- 
quark models [g[ 0| and for which the value of the ra- 
tio: M^^/m = 65, is definitive of light-quarks [Q. In 
addition, DCSB is also manifest in the vacuum quark 
condensate 



bo 



In 



IGeV^ 

^^QCD 



(0.221 GeV)^ 
(28) 



where we have used Aqcd = 0.2 GeV. The condensate is 
calculated directly from its gauge invariant definition |4q] 
after making allowance for the fact that Eqs. (p3|), (|2J) 
yield a chiral-limit quark mass function with anomalous 
dimension 7,„ = 1. This omission of the additional 
ln(p^/AQ(^j3)-suppression that is characteristic of QCD 
is a practical but not necessary simplification. 

Motivated by model DSE studies @, Eqs. (||), (|2|) 
express the dressed-quark propagator as an entire func- 
tion. Hence S{p) does not have a Lehmann representa- 
tion, which is a sufficient condition for confinement [Q. 
Employing an entire function for S(j}), whose form is 
only constrained via the calculation of spacelike observ- 
ables, can lead to model artefacts when it is employed 
directly to calculate observables involving large timelike 
momenta p9| . An improved parametrisation is therefore 



being sought. Nevertheless, no problems are encountered 
for moderate timelike momenta (see, e.g., Ref. [Q) and 
on the subdomain of the complex plane explored in the 
present calculation the integral support provided by an 
equally efficacious alternative cannot differ significantly 
from that of our parametrisation. 



2. Diquark Bethe-Salpeter amplitudes 

The renormalisation-group-improved rainbow-ladder 
DSE truncation, employed in Refs. ||, |l|, will 
yield asymptotic diquark states in the strong interaction 
spectrum. Such states are not observed and their ap- 
pearance is an artefact of the truncation. Higher order 
terms in the quark-quark scattering kernel (crossed-box 
and vertex corrections), whose analogue in the quark- 
antiquark channel do not much affect the properties of 
most of the colour-singlet mesons, act to ensure that 
QCD's quark-quark scattering matrix does not exhibit 
singularities that correspond to asymptotic (unconfined) 



diquark bound states |51 . Nevertheless, studies with 
kernels that do not produce diquark bound states, do 
support a physical interpretation of the masses obtained 



using the rainbow-ladder truncation, Eq. 



'99 



plays 



the role of a confined-quasiparticle mass in the sense that 
Iqq — l/nigq may be interpreted as a range over which 
the diquark correlation can propagate inside a baryon. 
These observations motivate the Ansatz for the quark- 
quark scattering matrix that is employed in deriving the 
Faddeev equation: 

[Mqq{k,q;K)Y;i = 



JP=0+A + , 



ik;-K)A^ {K)r^ {q-K). (29) 



While it is not necessary, one practical means of spec- 
ifying the T"^ in this equation, which is consistent with 
the above discussion, is to employ the solutions of the 
ladder- like quark-quark Bethe-Salpeter equation (BSE): 



{k;K) 



(2^ 



g{k-q)D'^l^{k-q) 



y^^lt.S{q + K/2)T' {q-K) 



^^,S{-q + Pl2) 



(30) 



where the effective coupling, Q (k) , is calculable using per- 
turbation theory for fc^ > 1 GeV^ and is modelled in the 
infrared (see, e.g., Refs. g |l|), and D^^^^'^ik) is the 
free gluon propagator. The amplitude is canonically nor- 
malised: 



2K^ = 



d 



tr 



d'^q 



dQ^ J (2^) 
X T{q-K)S^{~q + Q/2) 



f{q;-K)S{q + Q/2) 



(31) 
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Using the properties of the Gell-Mann matrices one 
finds easily from Eq. (|o|) that := F*^ C"^ satisfies 
exactly the same equation as the J colour-singlet meson 
hut for a halving of the coupling [Q. This makes clear 
that the interaction in the {qq)^^ channel is strong and 
attractive. The same analysis shows the interaction to 
be strong and repulsive in the {qq)e^ channel. 

A complete, consistent solution of Eq. ( ^0| ) requires a 
simultaneous solution of the quark-DSE, and while this 
combined procedure is not unmanageable it is a compu- 
tational challenge ^ In addition, we have al- 
ready chosen to simplify our calculations by parametris- 
ing Sip), and hence we follow Refs. ||, |ll |||5|] 
and also employ that expedient with F'^ , using the fol- 
lowing one-parameter forms: 

T"\k;K) - _i^i7«Cz75*r2.F(fcVwo^), (32) 
trl\k;K) = -l^i7-z7^Ct^.F(fcVO, (33) 

with the normalisation. A/'"' , fixed by Eq. (|3l|). Our 
Ansatze retain only that single Dirac-amplitude which 
would represent a point particle with the given quantum 
numbers in a local Lagrangian density: these amplitudes 
are usually dominant in a BSE solution ||2^, ^ |3^ , 



symmetry, we can focus on the A++, with it's simple 
flavour structure, because all the charge states are de- 
generate. The Dirac structure, though, remains complex 
and its general form is discussed in Rcf. Herein, as 
we have for the nuclcon, we use that study as the guide 
to a minimal model: 

-t+Fi^(ip[i2];i^)Ai:(i^)A.(^;P), (36) 

with 

A,(^; P) = S'^it, P) u,{P) ^+ + A^{£; P) ■ u{P) <f+ , 

(37) 

where Uu{P) is a Rarita-Schwinger spinor (see Ap- 
pendix ^, = £ -h P£ • P, and, again focusing on 
eigenfunctions of A+(P), 

5^(^;P) = ft{t\P)lD 

+ j^{il-i-i-PlB) f2ie;P), (38) 



A^{i;P) = [ati£;P)ID + ^atii;P)-^■i^]P^ 



(39) 



The Faddeev equation for the A now assumes the form 
Ap(fc;P) = 4 J ^M^Ak,e;P)A,{e;P) (40) 



3. Diquark propagators 

Solving for the quark-quark scattering matrix using the 
ladder-like kernel in Eq. (|30| ) yields free particle propa- 
gators for A'-' in Eq. (p9|). However, as already noted, 
higher order contributions remedy that defect, eliminat- 
ing asymptotic diquark states from the spectrum. It is 
apparent in Ref. |^ that the attendant modification of 
A'' can be modelled efficaciously by simple functions 
that are free-particle-like at spacelike momenta but pole- 
free on the timelike axis. Hence we employ |33| 



with 



AO^(if) 



1 



'0+ 



'1+ 



(34) 



'1+ 



where the two parameters mjp are diquark pseudopar- 
ticle masses and cujp are the widths characterising F'' . 
It is plain upon inspection that these Ansatze satisfy the 
constraints we have elucidated. 



B. Model for A 

The A is a spin-3/2, isospin-3/2 decuplet baryon and 
the general form of the Faddeev amplitude for such a 
system is complicated. However, as we assume isospin 



M 



— t^Fp (kq £qq/2]£qq)S (iqq kq) 
X t + f 1+ (iq - kqq/2; ~kqq) S{iq) (iqq) . (41) 



It is straightforward to construct four projection opera- 
tors that yield the coupled equations for /^2: 0,^2- 

We employ one more expedient to simplify our calcu- 
lations: we retain only the zeroth Chebyshev moments of 
/i,2, al^2: af^s; i-e-, we assume fi{£;P) = fi{P;P^), 
etc. We note that solving integral equations using a 
Chebyshev decomposition of the solution functions is 
a rapidly convergent scheme for isospin symmetric sys- 
tems ^ ^ and neglecting the other moments in 
this calculation will only have a small quantitative effect. 



C. Faddeev Equation Masses 

The nucleon and A masses can now be obtained by 
solving Eqs. (|o|), (^), and that also yields the bound- 
state amplitudes necessary for the calculation of the im- 
pulse approximation to A'^ and A form factors. The ker- 
nels of the equations are constructed from the dressed- 
quark propagator, and the diquark Bethe-Salpeter ampli- 



tudes and propagators, which are specified in Sees. [I A 1 



II A 3. These kernels involve four parameters. We fix 



mo+ = 0.74 GeV; 



(42) 
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TABLE L Calculated nucleon and A masses. The results 
in the first and third rows were obtained using scalar and 
pseudovector diquark correlations: mi+ = 0.90 GeV in row 
1, = 0.94 GeV in row 3. {mg+ = 0.74 GeV, always.) 

Pseudovector diquarks were omitted in the second and fourth 
rows, w/i 2 are discussed after Eq. (^^, and R in and after 
Eq. m8\). All dimensioned quantities are in GeV. 







Mn Ma 




0+ & 1+ 
0+ 


0.64 1.19 
0.64 


0.94 1.23 
1.59 


0.49 0.44 0.25 
0.39 0.41 1.28 


0+ & 1+ 
0+ 


0.45 1.36 
0.45 


1.14 1.33 
1.44 


0.44 0.36 0.54 
0.36 0.35 2.32 



i.e., we use the calculated scalar diquark mass in Eq. 
which is consistent with that obtained in recent, more 
sophisticated BSE studies (NB. TOo+ 2M^, 

Eq. (p7|), and hence it sets a good scale for nucleon 
observables.) This leaves mi+ and the diquark width 
parameters, ujjp . The immediate goal is to determine 
whether there are intuitively reasonable values of these 
parameters for which one obtains the nucleon and A 
masses: Mn = 0.94 GeV, Ma = 1.23 GeV, subject to 
the constraint mi+/mQ+ « 1.3, as in Eq. (||). 

The calculated masses are presented in Table |, from 
which it is apparent that the observed masses are eas- 
ily obtained using solely the dressed-quark and diquark 
degrees of freedom we have described above. The first 
two lines of the table also make plain that the additional 
quark exchange associated with the introduction of pseu- 
dovector correlations provides considerable attraction. In 
this case it reduces the nucleon's mass by 41%, in agree- 
ment with Ref. and, of course, without the l"*" corre- 
lation the A would not be bound in this approach. Fur- 
thermore, in agreement with intuition, the nucleon and 
A masses increase with increasing rujp . 

The values of the diquark width parameters are rea- 
sonable. For example, with 

Tp > Iq+ := i/i^o+ = 0.31 fm > Zi+ := I/0J1+ ~ 0.17 fm , 

(43) 

Tp is the proton's charge radius (experimentally, 0.87 fm), 
these correlations lie within the nucleon, a point also em- 
phasised by the scalar diquark's charge radius, calculated 
as described in Ref. M: 



rl+ = (0.55 fm)2 = (0.98 r^)^ 



(44) 



with calculated in the same model [Q. Further- 
more, defining w/^ 3 by requiring a least-squares fit of 
T{P' jiOj^ ^) to /i,2(^^), magnitude-matched at ~ 0, 
we obtain a scale characterising the quark-diquark sepa- 
ration 

:=lMi=0.40fm>0.15fm=^Zo+- (45) 



For the pseudovector analogue 



^9(99)ai = 0.36 fm > - 



(46) 



(a2(^^) is small in magnitude, slowly varying and not 
monotonic. Hence, in this case, the fit is of limited use. 
Nevertheless, it's momentum space width is roughly four 
times that of /i(£^).) 
For the A, 

^9(99) ,A = 0.35 fm K, = 0.32 fm > ^ Zi+ ; (47) 

and is important, characterised by a peak value of 
w -0.4 /f (£2 = 0) and w„a > 2tjyA, but is not: 

~ 0. 

The ratio 

R = /2 (^' = 0; -M2,)//i {P = 0; -M^^) (48) 

measures the importance of the lower component of the 
positive energy nucleon's spinor and it is not small, which 
emphasises the importance of treating these systems us- 
ing a Poincare covariant framework. For the A, R = 0.17. 

III. PION-INDUCED NUCLEON SELF-ENERGY 

We have illustrated that an internally consistent and 
accurate description of the nucleon and A masses is easily 
obtained using a Poincare covariant Faddeev equation 
based on confined diquarks and quarks. However, since 
the ttNN and irNA couplings are large, it is important to 
estimate the shift in the masses due to 7r-dressing. Herein 
we focus on the shift in the nucleon's mass because it is 
a much studied example. 



A. Model Field Theory: Linear Realisation of 
Chiral Symmetry 

We begin by considering a model n-N field theory de- 
scribed by the local Lagrangian density: 



C{x) = N{x)[i 

I J IT 



MV{x)] N{x) 



i^t,[d^V\x)d^V{x)\ 



16 



tr [2-V{x)-V\x)\ 



(49) 



(in this and Sec. Ill B we employ a Minkowski metric) 
where "tr" is a trace over Dirac and isospin indices, M 
is the nucleon's non-pion-dressed mass and the 7r-matrix 



1 



V{x) = exp ( i^^—f ■ Tr{x) 

J-rr 



(50) 
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with « 92MeV, the pion's weak decay constant. Ne- 
glecting the m^- (pion-mass-)term, this Lagrangian ex- 
hibits a hnear reahsation of chiral symmetry: 

N{x) N'{x) = V{ip)N{x), (51) 
N{x) N'{x) ^ N{x)V{ip) , (52) 
V{x) V{x) ^VHip)V{x)V\ip), (53) 

where V{ip) = exjp [i^r^r ■ ip / f^^) , with tp a spacetime- 
independent three- vector. (NB. The form of Eq. ( ^9| ) can 
be seen to arise from, and express, DCSB at the quark 
level using, e.g., the Global Colour Model ||. 
It also arises using the rainbow-ladder truncation of the 
DSEs.) 

Using Eq. (^) we explore the effect of 7r-dressing on 
M via the DSE for the nucleon self energy; i.e., S(P) in 

G-^{P) = P-M-J:{P). (54) 

In rainbow (Hartree-Fock) truncation that equation is: 



E(P) 



9 



^ A(A:^ 774) 75 G(P-fc) 75 



2M J (27r)4 
where, from the Lagrangian, 



9 



M 



(so that (7/1 = 1 at tree level in this model) and 
1 



fc2 



le 



1 



2 ujTlik) 



1 



1 



fco - uj^{k) 



fco + w^(fc) 



le 



(55) 



(56) 



(57) 



(58) 



with a;J(fc) = fc^ -|- m'^, is the free-pion propagator. The 



second contribution on the r.h.s. in Eq. ( p5[ ) is a tadpole 
(Hartree) term, which vanishes if the model is defined via 
dimensional regularisation. It is generated by the contact 
term in Eq. (p9|): /{2M) Nn ■ nN, whose presence and 
strength is dictated by chiral symmetry [ p5| , p6[ . 

As a first step we evaluate the self energy perturba- 
tively. To proceed with that we define the integrals 
in Eq. ( |55| ) by implementing a translationally invari- 
ant Pauli-Villars regularisation; i.e., we modify the tt- 
propagator: 



Aik\ml) ^ A^(fc2) = A(fc2, ) + ^ c, A(fc2, A^) , 

1=1,2 

(59) 

and then, with 

XI -ml 



Cl 



Xl ~ xl 



, C2 = 



A2 - Xl 



(60) 



Eq. 



yields 



A.(fc2) = Aik^ml) n (A^~m2)A(fc^A?), 



(61) 

in which case the integrals are convergent for any fixed 
A1.2. Furthermore, for TOtt ^ Ai ^ A2 = A 



A,(fc2) = A(fc2,m2)A2(fcVA2,l) 



(62) 



i.e., our Pauli-Villars regularisation is equivalent to em- 
ploying a monopole form factor at each ttNN vertex: 
g g A(fc2/A2, 1), where k is the pion's momentum [ p7[ . 
Since this procedure modifies the pion propagator it may 
be interpreted as expressing compositeness of the pion 
and regularising its off-shell contribution (a related ef- 
fect is identified in Refs. jH^, |5^) but that interpretation 
is not unique. 

In order to better understand the structure of the self- 
energy we decompose the bare nucleon propagator into a 
sum of positive and negative energy components 



G(P) 



G+{P) + G-{P) 



1 



ujn{P) 
+ A_(P) 



Pa-ujN[P)+ie 
1 



Po+ujN{P)-ie 



(63) 



(64) 



where uo%[P) = P^ + AP, and A±(P) (/'±M)/(2M), 
P = {uj{P),P), are, respectively, the Minkowski space 
positive and negative energy projection operators. Now 
the shift in the mass of a positive energy nucleon is 



:tr_D 



A+(P = 0)S(Fo = A/,P = 0) 



(65) 



We focus initially on the positive-energy nucleon's con- 
tribution to the loop integral; i.e., the A(fc)G"'"(P — k) 
contribution in the first term of Eq. (p5|), which we de- 
note by SpMl^. Evaluating the fco-integral by closing the 
contour in the lower half-plane, thereby encircling only 
the three positive-energy pion-like poles, we obtain 



SpMt 



-i9 

E 



■ (Pk u;jv(fc2)-M 



',=M.2 ^A,(fc2)K.(fc2)+C^Ar(fc2 



-M]' 

(66) 



with Co = 1, Ao = and w^. = fc + A^. It is obvious 
that SpM^ < 0; i.e., the Fock self-energy diagram's pos- 
itive energy nucleon piece reduces the mass of a positive 
energy nucleon. 

It is instructive to consider Eq. (^) further. Suppose 
that M is very much greater than the other scales then. 
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on the domain in which the integrand has significant sup- 
port, one has 



2M 



(67) 



and then 



SfMJ 



-3.9^ 



(27r)3 SAP ^ 

^ ' 1=0,1,2 



E 



so that 



(P 5fM+ 



r d?k P 
4M2 J (27r)3 

9 1 



1287r KP m. 
Thus, on the domain considered, 

3 9' 



(68) 

(69) 
(70) 



5fMX 



(71) 



where, as the derivation makes transparent, /(o,i) arc 
scheme-dependent functions of (only) the regularisation 
parameters but the first term is regularisation-scheme- 
independent. This first term is nonanalytic in the 
current-quark mass and its coefficient is fixed by chiral 
symmetry. (NB. If \i^2 are interpreted as setting a com- 
positeness scale for the ttNN vertex, and assume soft 
values; e.g. |§, - 600 MeV, then the quanti- 

tative value of SpM'^is completely determined by the 
regularisation-scheme-dependent terms.) 

We turn now to the A(fc) G^{P—k) contribution in the 
first term of Eq. (^5|), which we denote by SpM^. This 
describes the Z-diagram (anti-nucleon) contribution to 
the nucleon's mass and it is most eflticient in this case to 
close the fco-integration contour in the upper-half plane, 
thereby encircling only the three negative-energy pion- 
like poles: 



9 f d^ > 



(fk ujn{P)+M 
4u;jv(fc2) 



\ ^ Ci 

(72) 

It is obvious that SpMI^ > 0; i.e., the Fock dia- 
gram's anti-nucleon contribution to the positive-energy 
nucleon's mass is positive, and it is equally clear that, as 
evaluated with a pseudoscalar coupling Q , 



SfM+ + 5fM- = 3.9^ 



d^k M 



(2^)' 2ujN{k^) 



E 



^2 K.(fc2)+C.JV(A:2)]2-M2 

(73) 



i.e., dpM^ + SpM^ > 0, and hence that the negative- 
energy nucleon contribution overwhelms that of the 
positive-energy nucleon. HSfM+ > were the final word 
on the mass shift it would contradict all previous results 
for the effect of pion loops on the nucleon's mass ||6l| . 

Before addressing this issue we note that for M very 
much greater than the other scales then, reapplying the 
analysis that led to Eq. (|7l|), 



5fM^ = /(Q)(Ai,A2) + m^/(^j(Ai,A2) 
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327r2 M 



(Inm^ — 1) 



(74) 



The last term on the r.h.s. of this equation is an ad- 
ditional nonanalytic contribution to the nucleon's mass, 
and it is of lower order in 1 /M than the nonanalytic term 
in Eq. (0). This result, if it were to remain unamelio- 
rated, would also be in conflict with modern theory. 

Hitherto we have neglected the last term in Eq. (|5^), 
which describes the tadpole diagram's contribution to the 
positive-energy nucleon's mass shift, and the resolution 
of these apparent conflicts lies here. It is easy to evaluate 
and for Ai A2 = A 



ShAP 



3 9' 
327r2 AI 



[A^ + m2(lnK/AVl)] . (75) 



The inclusion of 6hAI^ solves both problems. It pro- 
vides for an exact, algebraic cancellation of the order- 
l/M term in SfAII^ that is nonanalytic in the current- 
quark mass, thereby ensuring that the nonanalytic term 
in Eq. ( [7l| ) provides the leading 0(1/Af) contribution to 
the nucleon's mass. The cancellation occurs because the 
(l/AI)^ contribution from the Z-diagram has the struc- 
ture of a tadpole term, for reasons that are intuitively 
obvious given that the (l/M)-expansion begins with an 
infinitely heavy nucleon. Furthermore, it must be exact 
because using dimensional regularisation, for example, all 
tadpole terms vanish and the leading nonanalytic term 
must be regularisation-scheme-independent. In addition, 
for aU A > 0, 



SAP = SfAII + SfAIV + ShAI+ < ; 



(76) 



i.e., the pion loop reduces the nucleon's mass. (NB. 
SAp{X) decreases monotonically from with increasing 
A, see Fig. 0.) 

We opened by asking for the scale of the mass shift 
produced by the pion loop. If we allow the interpreta- 
tion of the Pauli-Villars regularisation procedure as in- 
troducing a monopole form factor at each ttNN vertex, 
which modifies the pion's off-shell behaviour, then using 
soft values of the monopole scale: A ^ 0.5-0.7 GeV, as 
determined in quark-diquark Faddeev-amplitude models 
of the nucleon |^8[ ^ and inferred from data |60|, the 
0(g2) shift is as depicted in Fig. |l[ The magnitude is 
that of Refs. |, H. However, it is evident, and impor- 
tant to note, that this magnitude is extremely sensitive 
to the monopole's scale: centred on A — 0.6 GeV, a 10% 
change in A produces a 30% change in 6Ap. 
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0.50 0.55 0.60 0.65 0.70 0.0 0.5 1.0 1.5 

X (GeV) k (GeV) 



FIG. 1: Solid line: Shift in a positive-energy nucleon's mass 
due to the 0{g^) 7r-contribution to the self energy, Eq. ([r6|), 
obtained using a soft monopole pion-nucleon form factor to 
regularise the pion's off-shell behaviour. (M — 0.94 GeV, 
= 0.14 GeV and qa = 1.) SM+{X = 0.6 GeV) = 



-0.15 GeV. Dashed line: 



SpAI^; dot-dashed line: S' 



B. Nonlinear Realisation of Chiral Symmetry 

An alternative to Eq. ( ^9| ) is to build a Lagrangian 
density that contains only derivatives of the pseudoscalar 
field and thereby expresses a nonlinear realisation of chi- 
ral symmetry [ ]62[ . In chiral quark models such a La- 
grangian can be obtained via a unitary transformation of 
the fields in Eq. (B9)to obtain a so-called volume (pseu- 
dovector) coupling"! 54| . The leading term in the non- 
linear chiral Lagrangian can easily be obtained by using 
the equations of motion for a free nucleon to re-express 
Eq. ([49|). Neglecting that part of the Lagrangian density 
which describes the pseudoscalar field alone, this proce- 
dure yields: 



N{x) 



M -f-757'' r • 9p7f (x) -f . . . N{x), 



2M 



(77) 



and the rainbow truncation of this model's DSE is 
^2 I- d^k 



S(P) 



3i- 



(2^)4 



A(fc^mi)/^75G(P-fc)/^75. 

(78) 



No interaction survives that can generate a tadpole 
(Hartree) term. 

We again evaluate this self energy as a one-loop cor- 
rection to the positive-energy nucleon's mass. The con- 
tribution of the positive-energy nucleon is 



3.9^ 



d^k 1 



i=0,l,2 



16M2 J {2tt)3lun 
\j{ uJN - M) + 2p(u}x, + ^jv) 



, (79) 



FIG. 2: Illustrating that Pauli-Villars regularisation with 
finite mass-scales is a practical tool. Open circles - ii^(fc) 
calculated directly from Eq. ( ^l| ) using m.^ = 0.14 GeV, 
M = 0.94 GeV, = 1 and Ai ^ A2 = A = 0.6 GeV; solid 
line - least-squares fit to k^u^lP), which yields: u{k^) = 
0.99/(1 + PlX^f, A = 0.54 GeV; dashed line - limiting form 
from Eq. (p2). 



with cjjv = WAf(fc^), etc. Now, to make transparent the 
direct connection between our approach and other mass- 
shift calculations, we rewrite Eq. ( |79| ) in the form 



5^Af+ = 



67r 



fNN-n 



jfk Pu^{e) 

UJ^{P)[u;^{k^)+ujN{k^) - M] ' 

(80) 



where, as usual, /^jvir = g'^m^/{16TTM^) and, obviously. 



k'u\k') 



X J2 C ' 

1=0,1,2 



X^{lun - M) + 2k^{ujx^ + ujn) 



WAi +UJN - M] 



(81) 



This is useful because, for <C Ai — > A2 = A; i.e., 
on the domain in which Eq. (62) is valid, one finds alge- 
braically that 



u{k^) = 1/(1 + A:Va2) 



(82) 



which firmly establishes the qualitative equivalence be- 
tween Eq. ( |7^ ) and the calculation in Refs. p| p| 3^. 

In Fig. 1^ we compare the limiting form, Eq. (^), with 
u{k) calculated from Eq. (|l|). This emphasises the prac- 
tical utility of using a Pauli-Villars regularisation to rep- 
resent a ttNN vertex form factor. 

To provide a quantitative connection with other anal- 
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TABLE II: 5^M+ calculated using Eq. ^ in Eq. (gO|); 
i.e., a CBM estimate. The optimal bag radius for a descrip- 
tion of the neutron's electric form factor is J? = 0.95 fm = 
1/(0.21 GoV) 



R (fm) 


0.85 


0.95 


1.05 


-5'^M+ (GeV) 


0.091 


0.065 


0.048 



yses we employ 



u{k) 



M\k\R) 

\k\R 



(83) 



in Eq. (pO|); i.e., the CBM form for u{k), where R is the 
bag radius and ji {x) is a spherical Bessel function. The 
results are given in Table ^ and may be summarised as 
(in GeV) 



(0.065 ± 0.022) g\ , 



(84) 



where qa is the nucleoli's axial vector coupling constant. 
(NB. The result in Eq. (84) is also that obtained using 
a monopole form factor with the very soft scale Acbm = 
0.38 ±0.04 GeV.) We stress that in Eqs. (|9|) and (^) we 
used the couphng g = M/ /jr, Eq. (|56|), which corresponds 
to 5a = 1, whereas using the experimental value, gA — 
1.26, Eq. (H) gives 



'^CBM^-^I = -0.104 ± 0.035 GeV . 



(85) 



The larger shift described in Ref. |6| is obtained from 
Eq. (|8^) by using a smaller bag radius (~ 0.75 fm), which 
is needed to describe ttN scattering. The value of R em- 
ployed herein is appropriate to the calculation of nucleon 
electromagnetic form factors . A priori it is not clear 
which should be used for the calculation of hadron masses 



but recent lattice studies |27 favour a harder value. 

Returning to Eq. ( fz^ ) , the mass-shift contribution from 
the negative-energy nucleon; i.e., the Z-diagram, is 



2=0,1.2 



3g^ r (Pk 1 
16M2 J (27r)3 

XfjuJN + M) + 2fc^(^A. + ^at) 



(86) 



In this case we have S'^M+ < and 6^M- > but the 



sum 



3g2 f (fk M 



(87) 



8M2 ./ (27r)3 ujN 



E 



is self-evidently negative; i.e., with a pseudovector cou- 
pling the Z-diagram is much suppressed. 



Considering the heavy-nucleon limit again one obtains 
f 

327" 



^''mX = - 4^j^rni + /+.)(Ai, A2) mi -t- /+.)(Ai, A2) : 



(89) 

i.e., the same contribution, nonanalytic in the current- 
quark mass, as in Eq. (|7l|), but with different 
regularisation-dependent terms. In this case, however, 
because the Z-diagrams are suppressed by the pseudovec- 
tor coupling, the leading-order contribution to 5^M^ is 
0(l/A/)3. This is clear from Eq. (p6|), and makes imme- 
diately unambiguous the origin and nature of the leading- 
order nonanalytic contribution to the nuclcon's mass. 

Again interpreting the Pauli-Villars regularisation as 
introducing a monopole form factor at each ttNN ver- 
tex, we can estimate the magnitude of the 7r-loop's con- 
tribution to the nucleon's mass. Our results are depicted 
in Fig. 0. It is evident that S^M^ ^ SpM^, which 
illustrates the difference between the regularisation- 
dependent terms in Eqs. ( |7l| ) and (|8^). In addition, al- 
though it may not be immediately obvious. 



S^M+ = 6M- 



+ I 



(90) 



which is why there is only one solid curve in the fig- 
ure. This result provides a quantitative verification of the 
on-shell equivalence of the pseudoscalar and pseudovec- 
tor interactions, in perturbation theory, as long as the 
pseudoscalar interaction is treated in a manner consis- 
tent with chiral symmetry | ]56| . It also emphasises that, 
at least for estimating the mass shift, it is advantageous 
to employ the pseudovector interaction. We note, how- 
ever, that in fully embracing a Lagrangian density that 
expresses a nonlinear realisation of chiral symmetry one 
loses a direct correspondence with extant, ordered trun- 
cations of the DSEs, and hence also loses this correspon- 
dence between the Lagrangian's degrees of freedom and 
hadrons as composites of dressed-quarks. 



C. Model DSE 

We now build on the above analysis and seek a non- 
perturbative estimate of the 7r-loop's contribution to the 
nucleon's mass. Returning to the Euclidean metric de- 
scribed in Appendix ^ which is advantageous for numer- 
ical studies, the DSE for the nucleon's self energy using 
a pseudovector coupling is 



S(P) 



■ gly{P,k) ^^{{P ~ kf) 



(91) 



X7.(P-A:)75G(fc)7-(F-fc)75, 

with the following equivalent representations for the nu 
cleon propagator: 



G{k) 



l/[i"f ■ k + M 
l/[ij-kA(k^) 
—ij ■ k av{k'^) 



S(P)] , 

f M -f B{k^)] , 
^asik^), 



(92) 
(93) 
(94) 
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where M is the nucleon's bare mass, which is obtained, 
e.g., by solving the Faddeev equation. In Eq. (91), 
A7r(fc^) = l/[fc^ + m^] is the pion propagator, and 
gpv{P,k) is a form factor that we will use to describe 
the composite nature of both the pion and the nucleon. 
The self-consistent solution of Eq. ( |9l| ) yields A{k'^) and 
B{k'^), and thereby the nonperturbative mass shift. (For 
clarity we omit a discussion of renormalisation but re- 
mark on its effects following Eq. (119).) 
We now turn to the model specified by 



gpviP, k) = ^ exp(-(P - kf/A') . (95) 

The exponential form facilitates an algebraic evaluation 
of many necessary integrals and, as has been observed 
elsewhere [ p6| , is phenomenologically equivalent to a 
monopole form factor: 1/(1 + (P — fc)^/A^), if the mass- 
scales are related via A « V2A. Thus one can anticipate 
a quantitative correspondence between the A = 0.6 GeV 
monopole results of the previous subsections and those 
obtained in this with A « 0.9 GeV. 

Before proceeding with a nonperturbative solution of 
the nucleon's DSE we evaluate the one-loop self-energy 
so as to provide a direct Euclidean space comparison with 
Sees. Ill A , [II B . Using Eq. ( p5| ) we can evaluate the k^ 
integral to obtain 



A{t'') - 1 



3 a{t,^e-^^ M 

^ n I / — \ / — \ — 

327r^ M-^ Jq lUt^[k) ujn[k) 
327r2 M2 



(96) 



where a{t,K), b{t,K) are given in Eqs. (|aT9|)-(|A2i1). A 
and B are plotted in Figs. ||, |^. 

The one-loop-corrected nucleon mass: Af^ii , is the so- 
lution of 



M^.^2(- 



[M + B{- 



(98) 



and it is straightforward to show that Mjji — M = 
where (5Af+ is defined in Eq. (65). The calculated A- 
dependence of 5M^ is depicted in Fig. ||, and a com- 
parison with Fig. d reveals the equivalence between the 
Minkowski and Euclidean space formulations. 

The new feature in a nonperturbative study is that 
the position of the pole in the nucleon's propagator is 
not known a priori: locating it is the goal, and this pre- 
cludes an algebraic evaluation of the fc4-integral. The 
position of the pole will depend on the strength of the 
interaction and the nature of the form factor. In this case 
one must proceed by first evaluating the angular integrals 
in Eq. (|9l|), which are independent of G{k), noting that 
for a given function of {P ~ fc)^: 

J dilk f{{P - kf) = I y" dz^/l - z2 /(p2 + fc2 _ 2Pfcz) 

(99) 
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FIG. 3: Vector piece of the inverse dressed- nucleon propaga- 
tor. Dotted line: Ait) from Eq. (||), A{t = 4) = 1.001; solid 
line: numerical result for the one-loo p-dr ess ed fu nction in the 
spacelike region, obtained from Eqs. (IOC) - (103), which over- 
lies the dotted line in this region; das hed li ne: A it) obtained 
in the self-consistent solution of Eqs. (102), (103); dot-dashed 
line: Ait) {A{t — 4) — 1.00 2) obtained in t he self- consistent 
solution of Eqs. ( pdo| ), ( pdl| ) with Eqs. (Fl|), ([l|) added in 
the continuation to the timelike region. (All curves obtained 
with M = 0.94 GeV, = 0.14 GeV, = 1, A = 0.9 GeV.) 



This yields the kernels of the coupled, nonlinear integral 
equations for A, B: 



/C^(p2,fc2) 



/CB(p2,fc2) = 



dflkglviiP-kf) 



(P2 + fc2) 



{P^ - + ml{P^ + P) 



(P - fc)2 + ml 
2ml 



1 - 



(P - fc)2 + ml 



(100) 



(101) 



so that for spacelike p2 these integral equations can be 
written (a; = p2, y = P) 



x[A{x) -1] = - ^g^, 
Bix) = 



3 



dyylCA{x,y)crv{y), 



(102) 



dyyJCB{x,y) crs{y), 







(103) 



and solved numerically by iteration. 

To illustrate the accuracy attainable with thi s p roce - 
dure we: evaluated the integrals in Eqs. (IOC), (101) 
numerically for spacelike p2; inse rted Ajk'^) = 1 and 
S(fc2) = M on the r.h.s. of Eqs. (|lO|) and ( |l03| ); and 
calculated the integral over y numerically. This yields 
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FIG. 4: Scalar piece of the inverse dressed-nucleon prop- 
agator. Dotted line: B(t) from Eq. (||), M + B{t = 
4) = 0.937 GeV; solid line: numerical result for the one- 
loop-dressed function in the spacelike region, obtained from 
Eqs. (hOCl) - (103), which overlies the dotted line in this re- 



gion; das hed lin e: B {t) obtained in the self-consistent solution 
of Eqs. (Fo|), ([o|); dot-dashed line: obtained in the 

s elf-c onsistent solution of Eqs. (IIC), (111) with Eqs. (112), 
(113) added in the continuation to the timelike region. (All 
curves obtained with M = 0.94 GeV, = 0.14 GeV, gA = 1, 
A = 0.9 GeV.) 



the estimate of the one-loop-corrected nucleon propaga- 
tor in the spacelike region depicted in Figs. ||, ^. The 
agreement with the algebraic result is exact. 



The self-consistent solution of Eqs. (102), (103) in the 
spacelike region is easily obtained by iteration: the one- 
loop corrected functions are inserted on the r.h.s. to ob- 
tain the second iterate, which is then inserted on the 
r.h.s. to obtain the third iterate, etc., with the proce- 
dure repeated until the input and output agree within a 
specified tolerance. That happens very quickly, with the 
fourth iterate from free nucleon seed functions (^ = 1, 
B = M) agreeing with the third iterate to better than 
10^"* %. Hence "three pions in the air" are sufficient to 
fully dress the nucleon. The functions obtained in this 
self-consistent solution are also plotted in Figs. ^, |[ only 
A{t}^ is noticeably modified cf. the one-loop result. 

To locate the mass pole in t he n o npcr turbatively 
dressed nucleon propagator, Eqs. (102), (103) must also 
be solved for timelike . That req uir es a n analytic con- 
tinuation of the kernels in Eqs. (IOC), (101). The primary 
nonanalytic feature in their integrands is the pion pole 
and in continuing to timeUke it is necessary to prop- 
erly incorporate its effect. That is difficult when the ker- 
nels are only known numerically and an expeditious al- 
ternative is to develop an algebraic approximation, which 
is the approach we adopt. 

It is apparent that both kernels can be considered as a 
sum of two terms. The first is proportional to the angular 
average of gpy ((P— A:)^), and using Eq. (p5|) that integral 



FIG. 5: Dashed line: nucleon's one-loop mass-shift, calcu- 
lated from A, B in Eqs. (||), (|9^ - 5M+(A = 0.9 GeV ^ 
V2A) = -0.13 GeV; dot-dashed line: one-loop m ass shi ft ob - 
tained using the approximate kernels in Eqs. (105), (109); 
solid line: mass shift obtained via the self-consistent solu- 
tion of the nucleon's DSE using these approximate kernels - 
(5M+(A = 0.9 GeV) = -0.14 GeV. The dotted line is 5^MX\ 
i.e., Eq. ( |79| ) calculated in our Euclidean model. (All curves 
obtained with M = 0.94 GeV, = 0.14 GeV, qa = 1.) 



can be evaluated exactly: 
fpy{P',k') := / dnkglv{{P-kf) 



9 



AlVP 



-2(P^ + fe^)/A^ 



2Pk 



li{APk/K^ 



(104) 



where Ii {x) is a modified Bessel function and P = 

k — ^/k^. The second term in both cases is proportional 

to 



LUg2{P^,k^) := / dn 



{P -kY+ml ' 



(105) 



which, in general, cannot be expressed as a finite sum of 
known functions. However, if gpy is regular a.t P = k 
and its analytic structure is not a key influence on the 
solution, then the approximation 

iOg2{P\k^) « 9lv{\P^-k^\) fd^k 



{P - k f 



9lv{\P'-k^\). 



62 



(106) 
(107) 



where a = P"^ + k"^ + ml, b ^ 2Pk, is a reliable tool 
As these preconditions are obviously satisfied in our ap- 
plication - the dominant physical effect in ttN physics 
is the pion pole and that appears at a mass-scale much 
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lower than those present in gpy - we pursue our analysis 
using the following algebraic approximations: 



1 



9pv{x,y) (x + y) 



-9pv{x,y) 



2mi 



52 



(108) 



(109) 



To illustrate their efficacy, in Figs. H, ^ we plot the self- 
consistent solutions of 



x[A{x) - 1] = - 
B{x) = - 



167r2 



dyylCA{x,y) crviv), 



(110) 



dyyK.B{x,y) crsiv)- 



(111) 



The error introduced by the approximation is never more 
than 1% and is only that large for A{t'^ — 0). 

We can now define the model's analytic continuation 
to the timelike region. The approximate kernels' primary 
nonanalyticity is a square-root branch point whose ap- 
pearance and location are tied to the simple pole in the 
pion propagator, and in continuing to < q it is nec- 
essary to include the discontinuity across the associated 
cut. That is accomplished by adding the following 
additional terms to the r.h.s. of Eqs. (|ll(]| ), (111), respec- 
tively: 



167r2 



167r2 



dyy9'^{x,y)AICA{x,y)crv{y), (112) 
) 

dyy9'^{x,y)AJCBix,y)as{y), (113) 



where 



A/C^(x,y) = -^^^ [ix-y)'+ml{x + y)] 
2 mi 



(114) 



AitB{x,y) 



ml ^i^ + y + ^lY^^ (115) 



xy 



and y = Xb ^ —{\/—x — TOtt)^ is the location of the 
branch point. (NB. These terms are present only when 
P^+ml < 0.) The self-consistent solutions of Eqs. (110)- 
(113) are depicted in Figs. ||, ^ and unsurprisingly there 
is little difference between the one-loop results and the 
self-consistent solution. 

In Fig. ^we compare the exact one- loop mass shift with 
that obtained numerically using the approximate kernels. 



The error is never more than 5% with the approximation 
always overestimating the magnitude of the shift. (It is 
noteworthy that a large part of the one-loop mass-shift 
is due to the vector self energy; e.g., with A — 0.9 GeV, 
((5Af+)one-ioop is 40 % Smaller if the vector self-energy is 
neglected.) 

The fully dressed nucleon mass, Mu, is obtained by 
solving 



(116) 



with the nonperturbative mass shift given by = 
Mo — M . Again, this definition is completely equiva- 
lent to Eq. (^5|) evaluated at Mp, with the self-consistent 
solution of the DSE. The A-dependence of the nonpertur- 
bative shift is also depicted in Fig. ^ and comparison with 
the numerical one-loop result shows that the additional 
pion dressing adds < 5 % to |(5A/+|. 

Thus far we have used our Euclidean model to quanti- 



tatively reproduce the perturbative results of Sees. Ill A 



[II B and thereby make transparent the equivalence of 
the Euclidean and Minkowski formulations. In addition 
we have shown that the one-loop mass shift is ~ 95% of 
the total. 

However, we have not yet considered an effect of 
nucleon compositeness. A covariant ttNN vertex 
function must depend on three independent variables: 
gpy — gpY{P^, fc2^ (p _ ^-)2'j g^j-^i^ hitherto we have 

neglected its dependence on P^, k^. (We have al- 
ready seen that gpy — gpy{{P — fc)^) corresponds to a 
Pauli-Villars regularisation of the pion propagator alone.) 
The calculation of form factors that describe interac- 
tions between composite objects; e.g., studies of the p- 
u) mass splitting [ p9|, |69[ and electromagnetic form fac- 
tors 39, 1?^, ^l], indicates that the ttNN vertex 
should also suppress the pion-nucleon coupling when the 
nucleons are off-shell. We conduct an initial, exploratory 
study of this efi^ect by considering the product Ansatz 



gpy{P\k\P-k). 

_3_g-(P-fc)VA^ 

2M 



-(P^+M^+k^+M^)/Kl 



(117) 



which reduces to Eq. (^5|) when Ajv — > oo and guarantees 
gpy{-M'^,-M^,0) = g/{2M), as required. Previous 
applications of such a form factor in the ttN sector 
typically require 



Aat/A- 1.5 - 2.0. 



(118) 



(NB. While Aat is calculable using a covariant model of 
the nucleon, no such calculations exist and to constrain 
its value we must currently rely on phenomenology.) The 
effect on the mass shift of this off-shell suppression is 
depicted in Fig. ^: it is significant, leading to a reduction 
of > 50% in |(5M_|_|. For Ajv oo; i.e., in the absence of 
the off-shell suppression, this effect can be mimicked by 
a reduction in A; e.g., A ^ A' = 0.7 GeV yields 6M+ = 
-0.07 GeV, and we note that A'/v^ = 0.5 GeV, which is 
commensurate with Acbm ~ 0.4 GeV, after Eq. (84). 
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A, (GeV) 

FIG. 6: Effect on the mass-shift produced by incl udin g 
nucleon-off-shell suppression in the ttNN vertex: Eq. (117), 
for A = 0.9 GeV, = 0.14 GeV. (NB. SM+{A = 

0.9GeV,Ajv = oo) = -0.14 GeV.) 



Combining all the elements of our analysis we arrive 
at a result for the shift in the nucleon's mass owing to 
the TrA^-loop (for qa = 1.26, in GeV): 

-dM+ ~ (0.039 - 0.063) g\ = (0.061 - 0.099) . (119) 

In the preceding, for illustrative clarity, we did not 
account for the effects of finite vertex renormalisation; 
i.e., we set Zi = 1 = Z2 in Eq. (|9l|). Studies using 
the CBM indicate that a quantitative description of ttN 
vertex renormalisation requires that the A be treated 
on an equal footing with the nucleon and that this is 
crucial to obtaining a convergent expansion ||^, In- 
deed, one finds, as here, that the 7r-loop acts to sup- 
press the nucleon's wave-function renormalisation; i.e., it 
forces Z2 < 1, but in the CBM this effect is compen- 
sated by an almost matching suppression of Zi so that 
the bare and renormalised irN couplings are little differ- 
ent. A self-consistent, covariant treatment of the coupled 
composite- A^-A system is more than we are able to de- 
scribe herein. However, the CBM studies suggest that a 
reliable estimate of the effect of including the A can be 
obtained simply by solving an analogue of Eq. (^TJ) with 
Zi = Z2 for a renormalised model. 

We have done this and thereby arrive at a robust re- 
sult: the TrA^-loop reduces the nucleon's mass by ^10- 
20% Q. Extant calculations; e.g., Refs. ||, |, ||, show 
that the contribution from the analogous 7rA-loop is of 
the same sign and no greater in magnitude so that the 
likely total reduction is 20-40%. Based on these same 
calculations we anticipate that the A mass is also reduced 
by 77 loops but by a smaller amount (~ 50 - 100 MeV less). 

How does that affect the quark-diquark picture of 
baryons? To address this issue we again solved the 
Faddeev equations, this time requiring that the quark- 
diquark component yield higher masses for the N and 



A: Mat = 0.94 -I- 0.2 = 1.14 GeV, Ma = 1.232 + 0.1 = 
1.332 GeV. The results, presented in the third and fourth 
rows of Table ||, establish that the effects are not large. 
In this case omitting the axial-vector diquark yields 
AIn = 1.44 GeV, which signals a 10% increase in the 
importance of the scalar-diquark component of the nu- 
cleon. (It is an increase because this component now re- 
quires less correction. Note, too, that the scalar diquark's 
charge radius, ro+ = 0.63 fm, is 15% larger.) It also an- 
nounces a reduction in the role played by axial-vector 
diquark correlations in the nucleon, since now restoring 
them only reduces the nucleon's core mass by 21%, with tt 
self-energy corrections providing the remaining 14%. It is 
thus apparent that requiring an exact fit to the N and A 
masses using only quark and diquark degrees of freedom 
leads to an overestimate of the role played by axial- vector 
diquark correlations: it forces the 1+ diquark to mimic, 
in part, the effect of pions since they both act to reduce 
the mass cf. that of a quark-|-scalar-diquark baryon. 



IV. EPILOGUE 

We showed that an internally consistent description of 
the N and A masses is easily obtained using a Poincare 
covariant Faddeev equation that represents baryons as 
composites of a confined-quark and -diquark. We term 
this the "core mass" of the baryons. They are weakly 
bound in the limited sense that the sum of the masses 
of their primary constituents is little greater than their 
core mass. 

The on-shell irNN and TriVA couplings are large and 
hence it is conceivable that ttN and ttA self-energy cor- 
rections to the nucleon's mass may be significant. We 
therefore studied the effects of the TrA^-loop on the nu- 
cleon's core mass and found that, in well-constrained 
models, this loop reduces that mass by ^ 20%. Includ- 
ing the ttA self-energy contribution, the total reduction 
is likely to be between 20 and 40%. While this is a mate- 
rial effect it does not undermine the qualitative picture of 
baryons suggested by the Faddeev equation; namely, that 
baryons are primarily quark-diquark composites. This is 
consistent with the fact that a converged nonperturba- 
tive calculation of the 7r-induced self-energy requires only 
three "pions in the air," but to be certain we re-solved 
the Faddeev equation aiming at nucleon and A masses 
corrected for the tt self-energy contribution, and found 
little change in the character of the solution. 

The one notable effect was a material reduction in the 
nucleon's axial- vector diquark component. This is easily 
understood: ignoring 7r-loops forces the axial-vector di- 
quarks to mimic their effect. That surrogacy cannot be 
completely effective and may have led to quantitative er- 
rors, and errors of interpretation, in contemporary quark- 
diquark based calculations of quantities such as the neu- 
tron's charge form factor and the ratio /ipG^/G^j. Our 
results should serve as a signal of this possibility and 
stimulate increased caution and an objective reanalysis. 
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Our exploration of the role of 7r-loops was pedagogical. 
We made clear that the leading nonanalytic contribution 
to the miclcon's mass arises from that part of the loop 
integral which corresponds to a positive-energy nucleon; 
i.e., whether the ttA^'A'' coupling is pseudoscalar or pseu- 
dovcctor, the Z-diagrams do not affect the leading nonan- 
alytic behaviour. Furthermore, we showed explicitly that 
the one-loop mass shift calculated with a pseudoscalar 
coupling is precisely the same as that obtained with a 
pseudovector coupling, so long as, and only if, no dia- 
grams are overlooked in the pseudoscalar calculation. We 
illustrated that, using any translationally invariant reg- 
ularisation procedure which preserves information about 
the pion's finite size, the tadpole (Hartree) diagram gen- 
erated by a pseudoscalar coupling cannot be neglected 
because it balances the very large contribution from the 
pseudoscalar Z-diagram. This result should not be over- 
looked in the phenomenological application of model field 
theories founded on hadronic degrees of freedom. 
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APPENDIX A: EUCLIDEAN CONVENTIONS 

1. Metric and Spinors 
In our Euclidean formulation: 

4 



(Al) 



= 2 <5^^ ; 7^ = 7m ; '^/^^^ = ^ [7/^> 7^^] ; (A2) 
tr[757M7i'7p7CT] = -4 e^„pa- , £1234 = 1 • (A3) 
A positive energy spinor satisfies 
u{P, s){i-fP + M) = 0={ij-P + M) u{P, s) , (A4) 
where s = ± is the spin label. It is normalised: 

u{P, s) u{P, s) = 2M (A5) 
and may be expressed explicitly: 



i(P, s) = VM-i£ B-P 



M-i£ 



(A6) 



with £ = i^/WTAP, 



X+ ' / ' ^~ \ 1 



(A7) 



For the free-particle spinor, u{P, s) = u{P, 3)^74. 

Xhe spinor can be used to construct a positive energy 
projection operator: 



A4 



+ (^) ^= i E "(^' ^(^' ^) = ^ (-'^ -P + M). 



2M 



(A8) 



A negative energy spinor satisfies 

v{P, s) {i-i ■ P - M) ^0= {i-i ■ P - M) v{P, s) , (A9) 

and possesses properties and satisfies constraints ob- 
tained via obvious analogy with u{P, s). 

A charge-conjugated Bethe-Salpeter amplitude is ob- 
tained via 



r(fc;P) = CT{-k;PfC, 



(AlO) 



i=l 



where "X" denotes a transposing of all matrix indices and 
C = 7274 is the charge conjugation matrix, = —C. 

In describing the A resonance we employ a Rarita- 
Schwinger spinor to unambiguously represent a covariant 
spin-3/2 field. The positive energy spinor is defined by 
the following equations: 

(i7-P + M)M^(P;r) = 0, (All) 
l^u^{P-r) = 0, (A12) 
P^M^(P;r) = 0, (A13) 

where r = —3/2,-1/2, 1/2, 3/2. It is normalised: 

u^{P;r')u^{P;r) = 2M, (A14) 

and satisfies a completeness relation 

3/2 

— MP\r)u.{P;r)=K+{P)R^,, (A15) 

r=-3/2 
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where 



(A16) 



with — —1, which is very useful in simpUfying the 
positive energy A's Faddeev equation. 



2. Euclidean One-loop Calculations 

In Eqs. ^) and (|^ 

a(t, k) = -ao{t, n) + ai (t, n) + 02 [t, n) , (A17) 

6(i, k) = ao(i, k) +a2(t, , (AI8) 

where 

ao(i, k) = 



-Xo(f^w(t,K)) 



-^oi^Nit, k)) 



(A19) 



^oK(-))Ea:;r7;: 



f2Ar(i, K)Io{VLN{t, k)) 



(A20) 



a2{t,K) 



Q%it,n]-u;U^)y ^^^^^ 
where r^Ar(t, — cj_/v(/s:) + zt, $lA^(i, k) = ujn{i^) ~ and 



loix) := 2 / due-2"'/^'^^ 
= e2^'/'^'erfc(2;y2/A), 



A 

/2^ 



a;Xo(a:) 



(A22) 



(A23) 



where erfc(x) is the complementary error function and 
both these functions are odd under x ~x. 
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